Dynamic Fracture Model for Acoustic Emission 
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We study the acoustic emission produced by micro-cracks using a two-dimensional disordered 
lattice model of dynamic fracture, which allows to relate the acoustic response to the internal 
damage of the sample. We find that the distributions of acoustic energy bursts decays as a power 
law in agreement with experimental observations. The scaling exponents measured in the present 
dynamic model can related to those obtained in the quasi-static random fuse model. 

PACS numbers: 62.65.+k, 46.50.+a, 63.70.+h 



Crackling noise 0] is widely observed in systems as dif- 
ferent as superconductors Q, magnets || or plastically 
deforming crystals H . A typical example is the acoustic 
emission (AE) recorded in a stressed material before fail- 
ure. The noise is a consequence of micro-cracks forming 
and propagating in the material and should thus pro- 
vide an indirect measure of the damage accumulated in 
the system. For this reason, AE is often used as a non- 
destructive tool in material testing and evaluation. Be- 
side these practical applications, understanding the sta- 
tistical properties of crackling noise has become a chal- 
lenging theoretical problem. The distribution of crackle 
amplitudes follows a power law, suggesting an interpreta- 
tion in terms of critical phenomena and scaling theories. 
This behavior has been observed in several materials such 
as wood ||, cellular glass |J, concrete Q and paper || 
to name just a few. 

The statistical properties of fracture in disordered me- 
dia are captured qualitatively by lattice models, describ- 
ing the medium as a discrete set of elastic bonds with 
randomly distributed failure thresholds ||, O. Af- 
ter each crack the stress is redistributed in the lattice 
in a quasi-static approximation: i.e. the crack veloc- 
ity is much slower than stress relaxation. Thus acoustic 
waves are not taken into account and the activity is mon- 
itored by the damage evolution or by the dissipated elas- 
tic energy. Numerical simulations indicate that micro- 
cracks propagate in avalanches giving rise to an hetero- 
geneous response. The avalanche distribution is typically 
described by power law distributions and the results are 
usually interpreted in the framework of phase transitions 
0, |13|, 0, [16). Despite the fact that critical phe- 
nomena are normally associated with a certain degree of 
universality (i.e. the scaling exponents should not de- 
pend on micro-structural details), there has been so far 
no quantitative agreement between models and experi- 
ments. A reason that could account for this discrepancy 
is the absence of acoustic waves in most models. It is 
then not obvious how to relate AE activity to internal 
avalanches. 

Dynamic lattice models have been widely used in the 
past to analyze fracture processes [|l7| [l9|, , but al- 
though acoustic waves are explicitly included, the AE sig- 



nal is usually not analyzed. Here we use a lattice model 
for dynamic fracture in a disordered medium, to obtain 
a direct correspondence between the recorded AE activ- 
ity and the internal damage evolution. We find that the 
cumulative AE amplitudes are directly related — by a 
power law — to the cumulative damage. Next, we mea- 
sure the distribution of the AE burst energies and find a 
power law with an exponent j3 ~ 1.7 independent on the 
loading rate. This exponent can be related to the expo- 
nent describing failure avalanches in quasi-static models 

ii ii ii nm. 

We consider a scalar model of dynamic fractures where 
a two-dimensional lattice is loaded in mode III |18| : the 
lattice lies in the (x, y) plane and deformation occurs 
along the z axis, so that the equations of elasticity be- 
come scalar. The equation of the motion for the anti- 
planar displacement u of a site with coordinate i,j is 
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where the sum runs over the nearest neighbors (l,m) of 
site K is the elastic constant, p is the density and 

dissipation is simulated by a viscous damping with a con- 
stant r. In order to suppress some lattice effects, we use 
a 45 degree tilted square lattice. A constant strain rate 
is imposed to the model, by moving the boundary sites 
on two opposite boundaries at constant velocity V and 
— V, respectively. Periodic boundary conditions are im- 
posed in the other direction. Disorder is simulated as- 
signing randomly distributed failure threshold: a bond is 
removed (i.e. K is set to zero) when Am > / c , and f c is 
uniformly distributed in [0, 1]. Notice that in the quasi- 
static limit (V — > 0, p — ► 0, r — * 0) the model reduces to 
the random fuse model (RFM), where a lattice of fuses 
with random threshold are subject to an increasing volt- 
age |h| [ll| . Due to the scalar nature of our model there 
is a direct mapping between elastic and electric parame- 
ters p. 

The equation of motion (Eq. [I]) is integrated numer- 
ically using a fourth order Runge-Kutta method. We 
work with a lattice of linear size L = 80 and chose the 
units of space and time so that p — K = 1. Each time 
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FIG. 1: The stress-strain curve for different applied strain 
rates. 



a bond is stretched beyond its threshold the lattice con- 
stant is set to zero and an elastic wave is emitted. Due 
to the anti-plane constraint for the displacements, we 
only have transverse wave propagation with sound speed 
c = y/K/p = 1 in our units. The damping constant is 
chosen to be T = 0.1 so that typical length traveled by a 
wave is a little smaller than the lattice size. For smaller 
values of T ringing effects and reflected waves do not al- 
low to separate the single pulses and the lattice breaks 
at once. On the other hand, excessive damping leads to 
very small acoustic activity and the sample breaks sud- 
denly at the edges. Even if the damping constant is small 
reflected waves can induce boundary failure, due to the 
rigidity of the loaded edge. Thus we do not allow that 
bonds fail in two boundary layers of length I = 5 close 
to the loaded edges. This corresponds to apply a load 
through a soft contact. The model is simulated for a va- 
riety of loading velocities all much lower than the sound 
speed V <C c. 

Measuring the displacements of every lattice site and 
calculating the forces for every time steps, we have ob- 
tained the stress-strain curve for four different value of 
the applied strain rate. In Fig. |l| we show that the stress 
is a linear function of the strain up to the yield point, 
which precedes the total failure of the sample. The ap- 
plied strain rate has little effect on the linear part of the 
curve, while it influences the curve after the yield point. 

Monitoring the activity of some particular lattice sites 
we have direct access to the AE signal. These sites mimic 
the effect of transducers coupled to the material in a typ- 
ical AE experiment. In a typical run, we record the dis- 
placements, velocities and accelerations of four sites in 
the boundary layer and two sites in the interior. Typ- 
ically, AE distributions are recorded from a single site 
and averaged over ten realizations of the disorder. We 



have tested that the statistical properties of the signal 
do not vary for different boundary sites, while there is a 
clear difference between boundary and inner sites. In the 
following, we concentrate on sites in the boundary layer, 
in order to avoid excessive fluctuations due to failures 
occurring on neighboring bonds in the inner region. 

An example of the typical signals recorded with our 
model are reported in Fig. ||. A large acoustic activity is 
visible in the upper panel where we show the local accel- 
eration a of a boundary site as a function of time. We can 
also monitor the velocity signal which is simply related 
to the acceleration and display the same features. In the 
present model, it will be convenient to use the acceler- 
ation as a AE monitoring tool, since the velocity has a 
bias induced by the external loading: even in the absence 
of cracking the lattice has a non-vanishing velocity. We 
define the associated cumulative energy as 

E(t) = [ dt'a 2 (t'). (2) 
Jo 

The behavior of the cumulative acoustic energy E{t) is 
typically monitored in AE experiments. In some cases, 
E(t) is found to increase as a power law ||, or exponen- 
tially in other cases || . In general one expects a marked 
peak close to failure, as we also observe in Fig ^ ob- 
tained for V = 10~ 3 . The curve is well fitted by cubic 
law, E ~ t 3 . 
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FIG. 2: A typical signal measured in the model. The upper 
panel shows the velocity acceleration of a site close to the 
boundary. In the middle panel we show a magnification of 
a portion of the signal. In the lower panel is reported the 
cumulative energy, E(t) = J Q dt'a 2 (t'), as a function of time 
for a single realization of disorder. The dashed line follows t 3 . 

A central problem in AE measurements is to correlate 
the recorded acoustic activity with the internal damage 
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FIG. 3: The damage D evolution is displayed as a function 
of time for different loading rates. The curves collapse onto 
each other when plotted with respect to strain 7 = Vt/L. 



FIG. 5: The distribution of acoustic burst energies for dif- 
ferent driving rates. The line is a power law with exponent 
(3 = 1.7. 
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FIG. 4: The average cumulative energy plotted as a function 
of damage. The line represents a D 3 law which fits well the 
curve. 



state. In this way, AE can be used as a tool for damage 
evaluation. In our model, we have a direct access to 
the internal damage D that can be defined as the total 
number of failed bonds. We find that D increases linearly 
with time (see Fig. ||) apart from a rapid increase very 
close to failure. Rescaling the curves with the loading 
rate one sees that D is in fact a linear function of the 
applied strain 7 = (Vt)/L (see the inset of Fig. ||). 

These observations thus lead to a direct scaling relation 
between internal damage and released acoustic energy: 
Fig. [| shows that E scales as expected as D 3 . A direct 
consequence of this result is that the measured acoustic 



energy is proportional to the released elastic energy E e \ ~ 
KD~i 2 ~ D 3 ~ E. 

A large amount of theoretical activity has been devoted 
in the past to understand the origin of power law dis- 
tributions of AE amplitudes widely observed in material 
fracture. Most of the analysis was devoted to quasi-static 
models, such as the the RFM, where fracture was shown 
to occur in damage burst, distributed as P{D) ~ D~ T 
with r ~ 2.5 [jl2| [l3| . This value is in perfect agreement 
with the result r = 5/2 obtained exactly (2l[| for the ex- 
ponent of the avalanche distribution of the fiber bundle 
model (FBM) |22), where N fibers with random failure 
threshold are loaded in parallel. It was thus conjectured 
that the long-range stress transfer present in the RFM 
was equivalent to the infinite range load redistribution of 
the FBM, placing the two models into the same univer- 



sality class JL2], |13| . A similar exponent was found in a 
vectorial fracture model, so that this class could be even 



broader 13 . Comparing this result with AE experiments 
is problematic since quasi-static models do not account 
for wave propagation. 

Here we can directly measure the distribution of pulse 
sizes due to the acoustic activity. In Fig. || we report 
the distribution of energies for the acceleration signal, 
defining e = a 2 . In both cases the distribution decays as 
power law with an exponent (3 = 1.7±0.1, independent on 
the loading rate, which only affects the low part of the 
distribution. The same law is found in the case of the 
velocity signal. Experimental results report an exponent 
value in the same range, even if it differs a little from 
one material to another: for wood the exponent is (3 = 
1.51±0.05, for fiberglass (3 = 2.0±0.01 §, f3 = 1.30±0.1 
for paper B, {3 — 1.5 ± 0.1 for experiments on cellular 
glass @ . 

Using the scaling relation between released acoustic 
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energy and damage discussed above, we can relate the 
exponent j3 to r. From E ~ D 3 and D ~ t, we ex- 
pect e ~ -D 2 . Substituting this expression in the equa- 
tion for the probabilities P(e)de = P(D)dD, we obtain 
r = l+2(/3— 1) = 2.4, which is very close to t = 5/2 mea- 
sured in the RFM. Thus we conjecture that the acoustic 
energy exponent measured in our dynamic model is di- 
rectly related with the damage exponent measured in the 
corresponding quasi-static model |25j . 

In conclusions, we have introduced a lattice model of 
dynamic fracture which can be used to model AE exper- 
iments. The model allows to clarify important issues in 
the interpretation of the experiments, namely the rela- 
tion between internal damage and released acoustic en- 



ergy. In particular, we derive direct relations between the 
scaling behavior of failure avalanches and acoustic bursts. 
It would be interesting to generalize this analysis to more 
realistic situations, exploring the role of dimensionality, 
load conditions and lattice anisotropy. However, in com- 
paring the simulated signal with experiments, we should 
be careful about the definition of the events in the time 
series, since the amplifier and the AE sensors could bias 
the recorded waveform, introducing a systematic error in 
the data. 
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